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We describe an alternative approach to quantum computation that is ideally suited for today's 
sub-threshold-fidelity qubits, and which can be applied to a family of hardware models that includes 
superconducting qubits with tunable coupling. In this approach, the computation on an n-qubit pro- 
cessor is carried out in the n-dimensional single-excitation subspace (SES) of the full 2" -dimensional 
Hilbert space. Because any real Hamiltonian can be directly generated in the SES [E. J. Pritchett 
et al, arXiv: 1008.0701 , high-dimensional unitary operations can be carried out in a single step, 
bypassing the need to decompose into single- and two-qubit gates. Although technically nonscalable 
and unsuitable for applications (including Shor's) requiring enormous Hilbert spaces, this approach 
would make practical a first-generation quantum computer capable of achieving significant quantum 
speedup. 

PACS numbers: 03.67.Lx, 85.25.-j 
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I. INTRODUCTION 

In the standard gate-based model of quantum compu- 
tation with n qubits, the initial state of the (closed) sys- 
tem is represented by a 2"-dimensional complex vector 
if), and the computation is described by a unitary time- 
evolution operator U . Running the quantum computer 
implements the map 



ip — > Uip. 



(1) 



A key feature of this approach to quantum information 
processing is the exponential classical information stor- 
age capacity of the wave function ip. However, the num- 
ber of one- and two-qubit gates required to implement an 
arbitrary element of the unitary group U(2 n ) is at least 
2 2 ™, making the construction of an arbitrary U inefficient 
(using elementary gates). The goal of quantum algorithm 
design is to compute interesting cases of U with a poly- 
nomial number of elementary gates; two important ex- 
amples are the algorithms for factoring [l[ and quantum 
simulation [2.]. Although these algorithms are efficient, 
the number of gates required for interesting applications 
is still large and error-corrected qubits are therefore re- 
quired Q. 

Central to this gate (or quantum circuit) model of 
quantum computation is the idea that a set of elemen- 
tary one- and two-qubit gates can be used to construct 
an arbitrary high-dimensional U. But why should we 
use such a decomposition in the first place? The reason 
is that the Hamiltonians nature usually provides for us 
have one- and two-qubit terms (one-body and two-body 
operators). Consider, for example, the following some- 
what generic solid-state model of an array of n coupled 
qubits, 



H qc = e^c, + Y 9ij V °i 



i<3 



(2) 



written in the basis of uncoupled-qubit eigenstates 0). 
Here i, j = 1, 2, . . . , n, and fi, v = x, y, z. The e$ are the 
uncoupled qubit energies, gij (with i ^ j) are qubit-qubit 
interaction strengths, and is a fixed dimensionless 
tensor determined by the hardware. 

It is clear from ([2]) that one- and two-qubit operations 
are naturally generated in this system, and we know 
(from universality (a, l6j) that this is sufficient. If the 
hardware also provided controllable terms of the form 



ijk 



(3) 



then three-qubit gates could be used as primitives as well. 



II. QUANTUM COMPUTATION IN THE SES 

The idea we explore here is to perform a quantum com- 
putation in the rt-dimensional single-excitation subspace 
(SES) of the full 2 n -dimensional Hilbert space. This is 
the subspace spanned by the computational basis states 



|m) = cjjOO • • • 0) = |0- -0), 



(4) 



with m — 1,2,..., n. Experimentally, it is possible to 
prepare the quantum computer in the SES, and it will 
remain there with high probability as long as (i) the cou- 
pling strengths |<?y | are much smaller than the e*; and (ii) 
certain single-qubit operations such as 7r pulses are not 
used (however, 2-7T pulses are permitted and turn out to 
be extremely useful, and n pulses can be used to prepare 
SES states from the ground state 1 00 • • • 0)). 

The advantage of working in the SES can be under- 
stood from the following expression 0, H| for the SES 
matrix elements of model ((2]), 



2 



2J Z 



k<m 



9kr, 



9mk ) 

k > rn 



i<j 



jxx + jyy _ ^jxy _ jyxj 



(5) 



Therefore, we have a high degree of control over the part 
of the Hamiltonian in the SES. For example, in the simple 
case of an array of qubits coupled with a tunable erf <g> crj 
exchange interaction, 



quence, 



(m\H qc \m') = e m 5 mm i + g n 



(6) 



so the diagonal and off-diagonal elements are directly and 
independently controlled by the qubit energies and cou- 
pling strengths, respectively. Because of this high de- 
gree of controllability, n-dimensional unitary operations 
can be carried out in a single step, bypassing the need 
to decompose into elementary gates. This property also 
enables the direct quantum simulation of real but other- 
wise arbitrary time-dependent Hamiltonians, allowing a 
polynomial quantum speedup with sub-threshold-fidelity 
qubits 0. 

A quantum computer operating in the SES mode de- 
scribed here requires every qubit to be tunably coupled 
to the others, as implied by model ([2]). This requires 
n(n — l)/2 coupling wires and associated circuits. Con- 
sider, for example, the recently demonstrated tunable in- 
ductive coupler for superconducting phase qubits 0, [l(| 
(other tunable couplers for superconducting qubits have 
also been successfully demonstrated [TlTflij ). The circuit 
diagrams for a single phase qubit "q" and single coupler 
"c" are illustrated in Fig. [TJ where the crossed boxes rep- 
resent Josephson junctions. In terms of these elements, 
a possible layout for a fully connected array is shown in 
Fig. [5] for the case of n = 5. 

The SES computation method is not scalable, because 
it requires a qubit for every Hilbert space dimension used. 
However the advantage is that high-dimensional unitary 
operations can be carried out in a single step without 
the need to decompose the operation into elementary 
gates. This property allows processors of even modest 
sizes to perform quantum computations that would oth- 
erwise (using a gate-based approach) require thousands 
of elementary gates and therefore error-corrected qubits. 
In the remainder of this paper, we illustrate the SES 
computation method by applying it to Grover's search 
algorithm [15j ]. 



III. APPLICATION TO GROVER'S SEARCH 
ALGORITHM 

Here we consider the search for a single marked item 
in a database of size n, where n is the number of qubits 
in the SES processor ©. Within the SES, Grover's al- 
gorithm [15| is represented by the following control se- 



(WO m ,)" Grovcr t/unif|l) « \mf), 



(7) 



where t/ un ;f generates the uniform superposition of all 
SES states \m), m = 1,2, ... ,n. Here O m i is the oracle 
corresponding to the marked state \m'), W is Grover's 
inversion operator, and rtGrovcr ~ ( 7r /4)V« is the num- 
ber of iteration steps of the algorithm. In the following 
sections we show how each of these unitary operators can 
be generated in a single step. 



A. Preparation of the uniform state 

The uniform superposition state, 
|V\mif) 



v m— 1 v 



can be generated from the state |1) = [l 
where T stands for transposition, through 

|Vumf)=ie~e-^*-«|l), 



with 



and 



«unif = tt/ (2Vn) 



^unif = n/{2g\/n), 

using the SES Hamiltonian 

■2111. 
1 0. 
1 0. 
1 0. 



(8) 

■•] T . 
(9) 

(10) 
(11) 

(12) 



This can be seen by noticing that the spectrum of -ff un if 
and the transformation S (whose columns are shown un- 
normalized here for notational simplicity) that diagonal- 
s'^ -Hunif S are given by 



izes ff uni f via Hj^f 



and 
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FIG. 1: Definitions of two-terminal phase qubit circuit ele- 
ment "q" and coupler circuit element "c", expressed in terms 
of their conventional circuit diagrams involving inductors, ca- 
pacitors, current biases, and Josephson junctions (crossed 
boxes) [IIH. 




FIG. 2: (color online) Design for a fully connected network 
of n Josephson phase qubits. Each qubit "q" is suplied with 
a wire to which n — 1 coupler leads are attached. The case of 
7i = 5 is shown. The proposed design minimizes crossovers. 



B. Single-step oracle and W operators 

The oracle (corresponding to the marked state with 
ml = 3, for instance), is 







1 — 2|m )(m | 



10 
10 
0-10 
1 



(15) 



which can be generated via 



W = e~ lHot °, 



t = 7r/<5e, 



(16) 



using the SES Hamiltonian 



He 







-Se 





(17) 



8e is the detuning of qubit ml from the remaining qubits, 
all having common frequency e. This operation is simply 
a 27r rotation on the qubit associated with the marked 
state m! . (Although we have implemented it as a z rota- 
tion, an x or y rotation would work equally well.) Notice 
that in the limit g<E, which guarantees good isolation 
of the SES, the oracle can be made arbitrarily fast if we 
choose sufficiently large detuning St. 
Finally, the W operator, 

W = 2|VWf)Ww| - 1 

2-n 2 2 2 
2 2-n 2 2 
1 2 2 2-n 2 

// 2 2 2 2 -n 



(18) 



can be generated via 

% = (1 - n)ir/n, t w = ir/(ng), 



(19) 



with 



H w = g 



111 
10 11 
110 1 
1110 



(20) 



C. Comparison with gate-based computation 

respectively. Direct exponentiation then immediately Here we compare the SES and conventional gate-based 
leads to Eq. ©. The initial state |1) required in © approaches for a n = 256 item search. The search re- 
is easily prepared from the ground state via a 7r pulse. quires 12 iterations. 



4 



In the SES case, we choose a weak, g/2-K = 1.25 MHz 
coupling, which guarantees that the W operation is at 
least 1.5 ns long. We also choose a strong, 5c/2tt = 100 
MHz qubit detuning to guarantee that the oracle is suf- 
ficiently fast. Then t unl { — 12.5 ns, to — 5 ns, = 1-56 
ns, and the total duration of the algorithm is about 100 
ns, which is experimentally practical with current super- 
conducting architectures. 

The SES estimate should be contrasted with the con- 
ventional approach, which requires an input register of 8 
qubits, an output register of 1 qubit, plus 7 additional an- 
cilla qubits (to make the multiply controlled gates more 
efficient). The corresponding search oracle involves an 
8-fold CNOT gate, C 8 NOT, each of which can be made 
out of 85 two-qubit CNOTs (using 7 ancillas) Q. The W 
operator involves a 7-fold controlled-Z gate, C 7 Z, which 
uses 73 standard CNOTs (plus 6 ancillas). Thus, a con- 
ventional 256-item gate-based Grover search requires 158 
CNOT gates per search step. The full algorithm then 
contains nearly 2000 CNOT gates, plus local rotations. 

IV. DISCUSSION 

We have described an approach to superconducting 
quantum computation in which the computation is car- 
ried out in the single-excitation subspace of the full 



Hilbert space. Relative to the standard gate-based ap- 
proach, the SES method requires exponentially more 
qubits and is therefore nonscalable. The hardware re- 
quirements are also highly demanding: The fully con- 
nected array of n qubits requires n(n — l)/2 coupling 
wires and tuning circuits. But the SES approach is 
much more time efficient, permitting far larger computa- 
tions than currently possible using today's sub-theshold- 
fidelity qubits in the standard way. This was illustrated 
above for Grover's search algorithm. We would expect a 
similar result for Shor's algorithm, but note that practical 
factoring applications would require impossibly large pro- 
cessor sizes. General purpose time-dependent quantum 
simulation can also be carried out in the SES, allowing an 
n 3 polynomial quantum speedup Q- However, when a 
large scale error-corrected quantum computer eventually 
becomes available, the gate-based approach will perform 
better than the SES method. 
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